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Abstract 


The  Fourier  transform  method  is  applied 
to  the  problem  of  computing  viscous  flows 
involving  boundary-layer  separation,  based 
on  the  triple-deck  model  of  vlscous-lnvls- 
cid  flow  interaction.  As  used  here,  the 
metnod  is  pseudo-spectral  in  that  the  non¬ 
linear  inertia  terms  are  evaluated  in 
physical  variables,  although  the  main  com¬ 
putations  are  made  in  spectral  variables. 
The  Fast-Pourier-Transform  algorithm  is 
used  to  expedite  tne  Iterated  transforma¬ 
tions.  The  method  is  much  faster  than  con¬ 
ventional  finite-difference  methods;  typ¬ 
ically  only  ten  or  twenty  iterations  suf¬ 
fice  for  convergence  to  four  or  five 
digits  of  accuracy.  Furthermore,  no  arti¬ 
ficial  stabilization  schemes  are  necessary 
to  treat  the  reversed  flows  occurring  in 
separated  regions.  Results  are  presented 
for  both  incompressible  and  supersonic 
flows,  and  are  shown  to  compare  well  with 
previous  finite-difference  results.  > 


I.  Introduction 


In  recent  times  the  most  common  means 
of  solving  (numerically)  Prandtl'3  boun¬ 
dary-layer  equations  has  been  finite-dif¬ 
ference  marching  methods.  Despite  the  non¬ 
linearity  of  these  equations,  these  tech¬ 
niques  (incorporating  Newton  iteration) 
can  provide  rapid  accurate  solutions  on 
modern  computers .  The  marching  method  was 
applied  successfully,  for  example,  by  Jobe 
and  Burggraf1  in  their  solution  of  the 
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triple-deck  problem  of  interaction  between 
boundary  layer  and  free  stream  near  a 
trailing  edge.  Nevertheless,  there  are 
several  types  of  flow  for  which  these 
marching  techniques  either  fall  or  become 
difficult  to  implement.  Of  these  diffi¬ 
cult  types ,  perhaps  separated  flow  is  tne 
most  obvious.  The  physical  cause  of  the 
difficulty  is  the  region  of  reversed  flow, 
in  which  the  fluid  travels  in  the  direct¬ 
ion  opposite  to  that  of  the  main  body  of 
fluid.  As  a  result  there  is  a  change  of 
character  of  the  governing  boundary-layer 
equations.  Usually  these  are  of  parabolic 
type  for  which  marching  techniques  are 
valid;  however,  the  flow  reversal  due  to 
separation  changes  the  nature  of  the  boun¬ 
dary-layer  equations  to  quasl-elllptlc 
type,  with  information  being  propagated 
both  upstream  and  downstream.  Any  attempt 
to  obtain  an  accurate  solution  by  marching 
through  such  reversed-flow  regions  must 
fail  due  to  the  improper  direction  of 
Information  flow.  In  finite-difference 
procedures ,  the  failure  is  usually 
exhibited  either  by  lack  of  convergence  of 
the  iterative  procedure  or  by  severe  oscil¬ 
lations  in  the  "solution". 

Reyhner  and  Plugge-L5tz2  have  demon¬ 
strated  a  simple,  though  approximate, 
remedy  to  this  problem  of  treating 
reversed  flows.  Their  approximation  was 
to  neglect  the  product  of  the  streamwlse 
velocity  component  u  and  its  streamwlse 
derivative  ux  in  the  governing  equations 
wherever  u  became  negative.  This  techni¬ 
que  has  been  used  by  a  number  of 
3  b 

authors  ’  with  good  results  for  the  com¬ 
puted  pressure  distribution.  However, 
roucn  of  tha  success  of  this  technique  owes 
to  the  fact  that  the  reversed  flow  is  very 


slow  In  many  situations,  and  so  the  exclu¬ 
sion  of  this  component  of  the  Inertia  terms 
Is  of  little  consequence  In  those  cases. 

Williams-*  has  suggested  a  more  rational 
(and  more  complicated)  procedure  for  treat¬ 
ing  reversed  flows  Cased  on  bidirectional 
marching  techniques .  The  Reynner  and 
Fliigge-Lotz  (FLARE)  procedure  Is  used  as  a 
first  approximation,  followed  by  a  back¬ 
wards  sweep  In  the  reversed-flow  region 
only.  The  previously  neglected  uu^  term 
In  tne  reversed-flow  region  Is  now 
accounted  for  in  the  next  forward  sweep  by 
treating  it  as  a  known  quantity  evaluated 
from  the  preceding  backwards  sweep.  This 
downstream-upstream  Iteration  (DUIT)  Is 
then  continued  until  convergence  Is 
acnieved.  Williams  usually  finds  only  five 
to  ten  of  these  sweeps  to  be  necessary. 
Obviously  this  method  Is  considerably  les3 
straightforward  to  program  than  simple 
marching  procedures,  and  also  takes  rather 
longer  to  compute. 

Another  class  of  problems  for  which 
the  standard  marching  techniques  require 
Iterative  application  Is  tnat  to  free- 
lnteractlon  flows.  Originally  postulated 
by  Llghthlll**,  these  flows  are  essentially 
elgensolutlons  of  the  boundary-layer  equa¬ 
tions  that  provide  the  means  of  upstream 
Influence  In  the  flow,  despite  the  para¬ 
bolic  nature  of  the  governing  equations. 

7 

Stewartson  and  Williams'  and  later  Smith 

a 

and  Stewartson  used  a  shooting  approach 
to  generate  their  free-lnteractlon  solu¬ 
tions.  Here  the  pressure  Is  perturbed  by 
a  small  Jump  and  the  resulting  solution  is 
obtained  by  marching  downstream;  the  cor¬ 
rect  amount  of  perturbation  Is  determined 
such  that  the  free-lnteractlon  solution  Is 
suppressed  far  downstream.  As  usual  with 
shooting  methods,  tnis  approach  proves 
more  and  more  difficult  as  the  range  of 
Integration  is  extended  downstream,  owing 
to  the  ultimate  exponential  growth  of  the 
free-lnteractlon  solutions. 

The  above  methods  all  solve  the  steady- 
flow  equations  by  marching  In  the  flow 
direction.  An  alternative  Is  tlme-marchlng. 


with  which  the  unsteady  equations  are 
Integrated  forward  In  time  until  a  steady- 
state  solution  is  obtained.  This  approach 
permits  retention  of  all  the  inertia  terms 
In  the  equations  of  motion,  even  In 
reversed-flow  regions.  Jenson,  et  al^, 
and  later  Rizzetta,  et  al10,  applied  the 
time-marching  procedure  to  the  triple-deck 
formulation  for  supersonic  separated  flows 
over  ramp  configurations.  Their  method 
was  based  on  shear  stress  as  the  primary 
dependent  variable.  Upstream  Influence 
was  then  accounted  for  tnrough  an  Inter¬ 
action  condition  applied  as  a  wall  boun¬ 
dary  condl  ->n  on  the  stress:  namely,  the 
boundary-layer  compatibility  condition 
with  the  pressure  gradient  evaluated  from 
the  displacement  functlc  i.  This  condition 
Is  an  exact  requirement  of  tne  Interacting 
flow,  and  permits  downstream  conditions  to 
Influence  the  upstream  flow  at  each  time 
step  of  the  computation.  The  method  works 
well  but  flows  with  extensive  separation 
are  rather  expensive  to  compute  due  to 
their  slow  development  In  time.  An  Inter¬ 
acting  boundary-layer  program  based  on 
similar  Ideas  has  been  presented  by  Werle 
and  Vatsa^1  and  has  been  shown  to  give 
results  that  asymptote  the  triple-deck 
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results  for  very  large  Reynolds  number 

Various  alternatives  to  these  finite- 
difference  approaches  have  been  proposed. 
The  spectral  method  seemed  to  us  to  have 
characteristics  advantageous  for  computing 
separated  flows,  and  we  have  applied  It  to 
the  triple-deck  problem  in  this  study. 

The  (nonlinear)  governing  equations  are 
transformed  from  physical  to  spectral 
variables  using  the  Pourler  Integral  trans¬ 
form  In  the  main-flow  direction,  together 
with  finite-differences  In  the  transverse 
direction.  The  solution  is  computed 
(Iteratively)  in  spectral  space  and  then 
the  Inverse  transform  applied  to  obtain 
the  solution  in  physical  variables.  A 
major  advantage  of  the  method  Is  that 
reversed  flows  present  no  difficulty;  each 
point  In  spectral  space  relates  to  all 
points  on  the  path  of  Integration  In 


physical  space.  Thus  the  reversed  and  for¬ 
ward-flow  information  Is  diffused  together 
in  transform  varlaoies.  This  Interconnec¬ 
tion  between  the  physical  grid  points  via 
each  spectral  grid  point  endows  the  spec¬ 
tral  method  with  a  physically  implicit 
nature,  suggesting  an  accelerated  conver¬ 
gence  of  the  iterations.  As  will  be  seen, 
the  spectral  computations  exhibit  rapid 
convergence,  even  for  flows  exhibiting 
separated  regions  large  enough  that  other 
(finite-difference)  schemes  fall. 

Details  of  the  metnod  as  well  as 
results  of  computations  for  both  incompres¬ 
sible  and  supersonic  flows  are  given  below. 
In  both  cases,  the  spectral  results  are 
shown  to  compare  well  with  those  of  finite- 
difference  computations. 

II.  The  Mathematical  Problem 

The  problem  considered  Is  the  flow  dis¬ 
turbance  produced  by  a  shallow  perturbation 
of  the  surface  height  at  a  distance  L  from 
the  leading  edge  of  an  otherwise  plane 
wall.  Let  u,  be  the  flow  speed  and  the 
kinematic  viscosity  of  the  fluid,  the  sub¬ 
script  (-)  referring  to  conditions  In  the 
undisturbed  free3tream.  The  Reynolds  num¬ 
ber  R  «  u^L/v^  Is  regarded  as  large.  As 
Is  common  In  triple-deck  analyses.  It  Is 

convenient  to  define  the  small  parameter  e 

_  1/8 

as  R  .  Then  If  the  surface  height  per¬ 
turbation  Is  of  order  e3L  and  extends  over 
a  lengtn  of  order  e3L,  the  flow  disturbance 
Is  contained  in  the  triple-deck  structure 
originally  deduced  by  Messlter1^,  Nelland^ 
and  Stewartson^.  Since  this  asymptotic 
flow  structure  has  been  described  by  many 
authors,  we  merely  3tate  here  that  the 
problem  reduces  to  solving  the  conventional 
Incompressible  boundary-layer  equations 
(governing  the  lower-deck  flow)  subject  to 
an  unconventional  outer  boundary  condition 
whose  form  depends  on  whether  the  outer 
flow  Is  Incompressible  or  supersonic.  (A 
hypersonic  version  exists  as  well,  but  Is 
not  considered  In  this  work;  See  Ref.  15). 


It  is  convenient  to  describe  the 
triple-deck  disturbance  using  variables 
with  magnitude  of  order  one.  These  are 
denoted  by  upper  case  letters  while  lower 
case  letters  designate  the  corresponding 
physical  variable.  Thus,  the  coordinate 
representing  distance  parallel  to  the  sur¬ 
face  is  x  ■  e^aLX  while  the  normal  coor¬ 
dinate  Is  y  »  e^bLY.  The  origin  is  taken 
at  a  convenient  point  on  the  disturbed 
surface.  The  corresponding  velocity 
components  are  u  ■  c(d/b)u_U  and  v  « 
e^d/aJu^V,  while  the  pressure  p  »  p„(l  + 
e2cP).  Here  a,  b,  c  and  d  are  constant 
scale  factors  depending  on  Mach  number, 
wall-temperature  ratio,  and  the  surface 
stress  of  the  undisturbed  flow.  Their 
values  are  given,  for  example,  In  Ref3.  7 
and  10  for  supersonic  flow,  and  in  Refs.  1 
and  16  for  Incompressible  flow.  Following 

c  i  o 

Jenson,  el  al  ’  ,  we  Introduce  the  shear 

t  ■  as  basic  dependent  variable,  where 
subscripts  denote  partial  derivatives. 

Then  the  equation  to  be  solved  in  the 
lower  deck  Is  the  shear  transport  equation, 

UtX  +  VtY  “  TYY  (1) 

together  with  the  continuity  equation 

Ux  +  Vy  -  0  (2) 

The  conventional  no-sllp  conditions  apply: 

U  ■  V  ■  0  on  Y*0  (3) 

The  remaining  boundary  conditions  tatce  a 
fora  governed  by  the  triple-deck  structure. 
Upstream  the  shear  must  match  to  that  of 
the  undisturbed  flow  near  the  wall;  thus, 

t  -  1  for  X  -  --  (4) 

Moreover,  the  triple-deck  structure 
requires  that  the  main-deck  solution  cor¬ 
respond  to  a  simple  vertical  displacement 
of  the  original  undisturbed  boundary  layer 
by  the  lower  deck.  Consequently,  match¬ 
ing  of  main  and  lower  deck  solutions 
requires 

t-1  as  Y  ■*  “  (5) 


3 


The  corresponding  condition  for  U  is 
U  *  Y  *  A(X)  as  Y 

where  «(X)  is  cne  (negative)  displacement 
of  tne  main  deck.  For  supersonic  flow, 
tne  pressure  is  given  by  linearized  poten¬ 
tial  theory  as 

?  -  H 1  ( X )  -  A'U) 


t(u. 

Y) 

-  F(t (X, 

Y)) 

U(w , 

Y) 

»  F(U ( X , 

Y)} 

V(w, 

Y) 

•  F(V(X, 

Y)) 

In  terms  of  tnese  perturbation  variables, 
tne  shear  transport  equation  (1)  becomes 

Yrx  -  t y y  «  -(Utx  +  Vty)  3  R  (3) 


where  h'(X)  is  the  neignt  of  the  surface 
perturbation  In  lower-deck  scaling.  The 
ooundary-layer  compatibility  condition 
states  tnat  the  pressure  gradient  balances 
tne  normal  gradient  of  the  snear  at  the 
wall.  Expressing  A  in  the  form 


AU>  -  Lim  (U  -  Y) 

y— . 


(t  -  l)dY, 


the  compatibility  condition  for  supersonic 
flow  can  t.nen  oe  expressed  as 


t,.  -  H"  (X) 

1  Y-0 


rXX 


dY 


(5) 


On  the  otner  aand,  tne  elliptic  nature  of 
incompressible  flow  gives  rise  to  a  Hilbert 
Integral  for  tne  pressure: 

P(X)-i/"  H'(S)  ~  A-'(-S)  (7) 

3  -  X 


Thus,  an  interaction  condition  very  similar 
to  Eq.  (5)  can  be  given  for  incompressible 
flow.  Tnis  will  te  postponed,  however,  as 
the  similarity  is  mucn  more  remarkable  in 
spectral  variables. 

The  Fourier  transform  is  now  applied 
to  tne  governing  equations.  Using  an  over¬ 
tar  to  denote  transformed  variables,  tne 
transform  ♦(«)  of  any  real  function  b(X) 
is  given  as 

♦  (<-)  i  F( * ( X) )  3 


Actually  to  permit  convergence  of  the  inte¬ 
grals,  tne  transformation  i3  performed  on 
tne  perturoatlon  flow  quantities.  Define 

t  •  t  -  1 


where  the  second-order  perturbation  terms 
nave  been  placed  on  the  rignt  side  of  tne 
equation.  The  spectral  form  of  this  equa¬ 
tion  is 


IuiYt 


‘YY 


R  =  - (Ul ^  +  Vty) 


(9) 


The  rignt-side  function  R  has  the  form  of 
a  complicated  convolution  integral,  how¬ 
ever,  it  need  not  be  given  here,  since  the 
term  will  be  evaluated  by  another  method 
(pseudo-spectrally )  .  The  continuity  equa¬ 
tion  (2)  is  expressed  spectrally  as 

r  Y 


-lwU 


-lu 


tdY 


(13) 


The  boundary  conditions  ( 3 ) - ( 6  )  also 
must  be  expressed  in  spectral  variables. 
Thus, 


U  - 

0 

on 

Y  -  0 

(11) 

V  - 

0 

on 

Y-0 

(12) 

T  ♦ 

0 

as 

Y  ♦  - 

(13) 

and  the 

interaction 

condition  (6) 

for 

supersonic  flow  becomes 


-w^(H(w)  - 


Y-0 


t(w,  Y)dY)  (19) 


For  incompressible  flow,  the  pressure 
is  given  oy  Eq.  (7),  so  tnat  the  compati¬ 
bility  condition  becomes 


luF 


Y»0 


Applying  tne  Fourier  transform  to  tne 
Hilbert  Integral  in  Eq .  (7)  thus  yields 
tne  Incompressible  interaction  condition 


r- 


Y-0 


-lw  |  u  |  ( H  ( w  )  - 


t  ( u> ,  Y)dY)  (15) 
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TtMiii  i  .  _  AaMMiaaAfli - 1 


U  -  U 


Y 


wnlcn  Is  very  similar  In  form  to  the  super¬ 
sonic  Interaction  condition.  In  fact,  this 
slignt  difference  between  Eqs.  (14)  and 
(15)  Is  the  only  difference  In  the  triple- 
deck  theories  for  supersonic  and  Incom¬ 
pressible  flow.  Yet  as  will  be  demon¬ 
strated,  large  apparent  differences  In  the 
pnyslcal  flow  properties  result.  The 
significance  with  regard  to  computer  pro¬ 
gramming  Is  that  tne  change  of  a  single 
FORTRAN  statement  permits  computing  either 
supersonic  or  Incompressible  flows  with 
tne  same  computer  code. 


III.  The  Solution  Procedure 

The  3oectral  shear-transport  equation 
(9)  nas  been  written  In  the  form  of  a 
linear  ordinary  differential  operator  on 
tne  left  side,  with  the  non-linear  inertia 
terms  collected  together  Into  the  right- 


side  function  3. 


7he  solution  Is  obtained 


oy  Iteration,  solving  flr3t  with  R  set 
Identically  zero;  then  In  later  Iterations 
R  Is  evaluated  from  the  subsequent  solu¬ 
tion,  as 


-  ?(n) 
-  Tr{ 


"(n-1) 


-  -(D(n“l)ij‘n-l>  +  v(n-1)T^n"1)) 

Central  differences  In  Y  are  applied 
to  Eq.  (lc),  as  well  a3  to  the  interaction 
condition,  either  (14)  or  (15).  The 
result  of  tne  first  Iteration  Is  the 
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linearized  tneory  of  Stewartson  ,  which 
provides  a  reasonably  accurate  prediction 
for  small  disturbances.  For  disturbances 
large  enougn  to  produce  flow  separation, 
furtner  Iteration  Is  necessary  to  bring  in 
tne  effects  of  the  non-linear  Inertia 


The  right-side  function  R  In  Eq.  (16) 
can  ne  expre  ised  as  a  convolution  of  the 
velocity  referred  to  two  different  points 
In  spectral  space.  However,  evaluating 
tals  convolution  Integral  Is  a  very  Inef¬ 
ficient  way  of  determining  R,  requiring 
of  order  multiplications,  for  each 


complete  iteration,  where  .Jy  and  are 
the  number  of  grid  points  In  the  Y  and  u 
(spectral)  directions,  respectively. 

A  much  more  efficient  method  is  to  use 

tne  so-called  Fast  Fourier  Transform  (FFT) 

13 

algorithm  of  Cooley  and  Tukey  .  The  pro¬ 
cedure  used  here  Is  to  invert  tne  trans¬ 
forms  t  and  Iut  to  obtain  the  physical 
variables  t  and  tx<  From  these  U,  Ux>  and 
V  are  determined,  permitting  the  inertia 
function  R  of  Eq.  (8)  to  be  evaluated. 
Applying  tne  direct  transform  to  R  then 
results  in  the  desired  Inertia  transform 
tern  R.  Eq.  (16)  is  then  solved  for  tne 
next  approximation,  and  the  above  procedure 
Is  iterated  until  tne  solution  repeats  to 
tne  desired  number  of  decimal  places.  The 
advantage  of  this  method  arises  from  the 
efficiency  of  the  FFT  algorithm;  tne  scheme 
outlined  here  requires  only  of  the  order  of 
log2  multiplications  for  each  com¬ 
plete  Iteration,  compared  with  Hy.'J2  for  tne 
direct  metnod  based  on  the  convolution  Inte¬ 
gral.  An  order  of  magnitude  reduction  In 
computing  tine  Is  achieved  for  quite  rea¬ 
sonable  grid  resolution.  It  may  be  noted 
that  tne  method  used  here  is  the  reverse 
of  the  pseudo-spectral  metnod  described  by 
Orszag1^  and  Roache20,  wno  solve  tne  equa¬ 
tions  of  motion  In  the  physical  plane  ana 
use  tne  Fourier  transform  to  evaluate 
derivatives . 

The  discretization  of  the  Fourier 

transforms  was  carried  out  as  follows. 

Define  the  discrete  physical  and  spectral 

variables  X,  and  w.  as 

J  * 


Xj  -  (J-l-M)dX  for  J  «  1,2,...,:;;  (17a) 


«k  -  (k-l-M)Au  for  k  -  1,2,...,:;;  (17b) 


where  for  convenience  we  take  N  •  2M. 
Denote  any  physical  function  d(Xj)  oy  «j, 
and  tne  corresponding  spectral  function 
$(wk)  by  5^.  Choose  the  range  of  vari¬ 
ables  sucn  that  »(X)  is  negligable  for 
X  <  X^,  X  >  and  5(<w)  Is  negligable 


for  w  <  ,  u>  > 


Then  tne  Fourier 
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Integral  transform  b(w),  defined  In  the 
last  section,  can  be  approximated  by  the 
finite  sum 

Ji  -1M. 

e  J  (18) 

J-l 

The  range  of  Integration  Is  slightly 
uncentered,  from  X  ■  -(M  +  1/2)AX  to 
X  ■  (M  -  1/2)AX.  Since  b(X)  Is  real,  it 

follows  from  the  form  of  the  transform 

-  -« 

tnat  b(-w)  «  «  (w),  where  the  asterisk 
denotes  the  complex  conjugate.  3ecause  of 
this  property  the  transform  variables  need 
not  be  stored  for,  say,  positive  values  of 
The  midrange  parameter  M  In  Eqs.  (17) 
Is  now  required  to  be  an  even  integer  and 
the  grid  spaclngs  AX  and  Au  satisfy  the 
relation 


Then,  Eq.  (13)  can  be  reduced  to  the 
finite  Fourier  transform 

JL  ,  -12s (k-l) (J-l)/N 

♦k  ■  e  (20) 

J-l 

For  the  Inverse  Pourier  Integral  trans¬ 
form,  the  centered  range  of  Integration 
•  -MAui  to  -  MAu  is  taken  and, 

noting  that  b,j+1  -  the  trapezoidal 
rule  results  In  the  finite  sum 


*j  •  §7  ci(^i  +  ♦*>  ♦  21  *k  e  k 


*J  -&[(-  1)J  §££>1)^ 


12s(k-l ) ( J-l )/U 


Au  rl,"» 
7?  l-7(*i 


5k  e  k  J3 


Since  Is  real  and  (b^  -  b^)  Is  Imagin¬ 
ary,  the  result  Is 

♦j  t  5k  ^ 

k-1 

Substituting  for  and  Xj  from  Eqs.  (17) 
and  (18),  and  again  recalling  that  M  Is 
even,  the  Inverse  finite  Pourier  trans¬ 
form  results: 


Eqs.  (20)  and  (21)  are  best  evaluated  using 
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the  Cooley-Tukey  FFT  algorithm  ,  for  the 
reasons  discussed  above.  In  that  case,  A 
must  be  a  power  of  2. 

The  range  of  integration  was  approxi¬ 
mately  centered  about  the  origin  In  both 

physical  and  spectral  variables  because  of 
—  -# 

the  symmetry  property  »(-w)  »  b  (w).  The 
real  variable  b(X)  may  be  quite  unsymmet- 
rlcal,  however,  and  aliasing  errors  may  be 
observed.  This  effect  arises  because  the 
non-periodic  function  *(X)  Is  represented 
by  the  finite  Fourier  transform  as  the 
periodic  function  b(X),  with  period  MAX. 
Consequently,  If  an  Improperly  aligned 
interval  of  length  .VAX  is  sampled  from  b(X), 
aliasing  may  occur  on  either  the  left  or 
right  side  of  the  sample  interval.  If  on 
the  left,  for  example,  the  aliased  values 
actually  would  correspond  to  the  correct 
values  taken  from  b(X)  on  the  right  of  the 
Interval.  No  such  errors  exist  In  b(u)  as 
computed  here,  since  the  Inertia  terms  are 
evaluated  at  each  X-grid  point  Indepen¬ 
dently.  Aliasing  does  occur  In  the 
physical  results,  but  Is  easily  recognized 
since  the  flow  properties  decay  to  the 
undisturbed  31aslus  values  at  infinity. 

The  effect  can  be  avoided  either  by 
restricting  the  X-range  of  the  data  pre¬ 
sented,  or  by  use  of  a  special  Inversion 
formula  based  on  a  properly  uncentered 
Integration  Interval.  Both  metnod3  were 
used  for  the  results  presented  In  the  fol¬ 
lowing  section.  A  better  method  of  sup¬ 
pressing  aliasing  errors  would  be  to 
evaluate  the  truncated  tails  of  the  Fourier 
Integral  by  asymptotic  analysis  and  add 
their  contributions  to  the  FFT  result. 

This  procedure  Is  cumbersome  and  was  not 
necessary  for  the  cases  dealt  with  here. 

IV.  Results 

Two  surface-height  configurations, 
Indicated  by  H^(X)  and  H2(X),  have  been 


chosen  to  Illustrate  the  spectral  triple- 
deck  metnod.  The  shapes  are  shown  In  Fig¬ 
ure  1,  ana  are  defined  by  the  relations 

H1(X)  -  a/(l  +  X2) 

!  0  for  |  X  |  >1 

h2(X)  -] 

(o(l  -  X2)2  for  |X|  <  1 

The  first  of  these  snapes  exnlblts  a  long 
gentle  variation  of  height,  but  has  a 
Fourier  transform  that  Is  much  sharper, 
being  significant  only  over  a  narrow  range 
of  ui : 

»  cnre”^ 

The  second  shape  has  the  opposite  behavior: 
it  differs  from  zero  only  over  a  narrow 
range,  but  has  a  slowly  varying  Fourier 
transform: 

-  2  5 

H2(u))  «  -l6[3i»>cos  u)  +  (w  -  3)  slnu]AiJ 
wltn 

H2(0)  -  16/15 

The  opposing  behavior  of  these  two  shapes 
In  physical  and  spectral  coordinates  make 
them  useful  for  case  studies  of  the  spec¬ 
tral  method.  The  second  snape  H„(X)  was 
21  ‘ 

used  earlier  by  Sykes  ,  and  his  finite- 
difference  results  will  serve  as  a  check 
case  for  Incompressible  flow.  For  super¬ 
sonic  flow,  an  Independent  check  case  was 

run  for  H.(X)  using  the  tlme-marchlng 
1  9 

finite-difference  program  of  Jenson  and 

Rlzzetta^ ,  modified  to  Incorporate 
3econd-order  accurate  Crank-Nlcolson  dif¬ 
ferencing  In  X. 

Figure  2  displays  results  for  super¬ 
sonic  flow  over  tne  surface  contour  H^X) 
for  values  of  amplitude  a  from  -3  to  ■*■3. 
(Positive  a  corresponds  to  a  hump,  nega¬ 
tive  a  to  a  hollow  In  the  surface.)  The 
curves  represent  surface  stress  computed 
by  the  spectral  program  for  *  128, 

AX  -  0.25,  Ny  -  25,  AY  ■  0.5.  The  grid 
Interval  Is  obvious  in  the  figure  since 
values  at  grid  points  are  connected  by 
straight  lines.  The  symbols  are  for  the 


checx-case  results  from  tne  finite-dif¬ 
ference  program  for  a  »  -  2,  with  .'1^  ■  81, 
AX  -  0.25,  My  -  37,  AX  -  0.25.  Obviously, 
the  results  of  the  two  metnods  are  in 
excellent  agreement.  The  spectral  method 
was  much  faster,  however,  requiring  only 
10  iterations  for  a  *  -2,  while  the  time¬ 
marching  finite-difference  program  required 
300  time  steps  to  converge  to  the  steady 
state.  In  terms  of  processing  time,  tne 
spectral  method  Is  an  order  of  magnitude 
faster,  requiring  for  these  cases  about  10 
seconds  CPU  time  on  an  Amdahl  A70/V6  com¬ 
puter  verses  about  81  seconds  for  tne 
finite-difference  program. 

For  the  humps  (a  »  1,2,3),  the  flow 
coming  from  the  left  at  first  decelerates , 
reducing  the  shear  below  the  Blaslus  value 
(t  »  1  In  lower-deck  scaling).  For  suf¬ 
ficiently  tall  humps  (a  >  7,  approx.)  the 
flow  would  be  expected  to  separate  here, 
as  well  as  on  the  lee  3lde  of  the  hump. 

Near  the  top  of  tne  hump  the  lower-deck 
flow  is  squeezed  by  the  rapidly  Increasing 
surface  height  and  the  surface  stre33 
peaks  to  very  high  values.  For  the  hol¬ 
lows  (a  *  -1,-2, -3)  the  reverse  of  these 
effects  occurs.  In  this  case,  a  flow- 
separation  bubble  occurs  at  the  bottom  of 
the  hollow  for  a  <  -2.6  (approx.).  The 
non-linear  effects  can  be  seen  In  Figure 
2,  but  are  more  clear  In  the  table  of 
values  below,  where  the  difference  between 
tne  center-point  shear  and  the  Blaslus 
value  of  unity  Is  normalized  with  respect 
to  a.  The  effect  of  non-linearity  Is 
demonstrated  by  the  deviation  from  the 
value  0.537  for  a  ■  0.  Clearly  non¬ 
linearity  is  more  important  for  the  hol¬ 
lows  than  for  the  humps,  accounting  for  a 
change  of  about  30*  from  the  linear-theory 
value  for  a  •  -3. 

The  effect  of  spatial  resolution  on 
accuracy  Is  indicated  In  Figure  3  for 
a  *  -2.  The  symbols  represent  results 
for  a  coarse  grid  having  N  -  32,  AX  -  1.0, 
and  the  curve  represents  fine-grid  results 
with  Nu  -  128,  AX  •  0.25.  (The  fine-grid 
results  are  the  same  as  those  of  Figure  2.) 
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For  both  grids  !ly  -  25  and  AY  -  0.5.  Both 
cases  used  the  same  X-range  (-16  <  X  S  16 , 
approx.),  whlcn  Is  clearly  large  enough  to 
yield  accurate  results.  As  N  is  reduced 

OJ 

below  32,  with  the  same  X-range,  the 
accuracy  degrades  quickly  owing  to  poor 
approximation  to  tne  surface  shape  for 
AX  >  1.  Also  reducing  tne  X-range  below 
about  -6  X  <_  3  severely  degrades  the 
accuracy,  most  likely  due  to  truncation  of 
tne  wake-lixe  flow  decay  downstream. 


TA3LE  I 

SURFACE  STRESS  AT  X  *  0 
H<  X)  -  a/(l  +  X2) 


A  comparison  of  separated  supersonic- 
flow  results  with  those  for  incompressible 
flow  past  the  surface  contour  H^(X)  is 
given  in  Figure  4.  The  amplitude  a  •  -5 
was  cnosen  for  the  comparison,  since  the 
Incompressible  flow  does  not  separate  for 
a  <  -4.  The  incompressible  flow  (upper 
curve)  converged  in  30  Iterations,  but 
the  supersonic  flow  (lower  curve)  did  not 
converge  to  the  same  accuracy  until  80 
iterations.  In  the  supersonic  case,  the 
residuals  did  not  dec^y  monotonically , 
and  an  under-relaxation  factor  of  0.5  was 
applied  to  obtain  convergence.  This  case 
wa3  the  only  one  of  those  analyzed  in  thi3 
3tudy  that  required  under-relaxation.  In 
general,  the  supersonic  and  incompressible 
flows  show  quite  similar  surface  shear- 
stress  patterns,  especially  in  the  sepa¬ 
rated  region.  It  may  be  noted  that  the 
separated  t0(X)  distribution  13  very  much 
like  that  computed  for  other  configura¬ 


tions  such  as  In  the  ramp  study  of 
Rlzzetta,  et  ali0. 

Our  last  comparison  is  for  tne  quartic 
hump  H p ( X )  ,  which  provides  a  severe  test 
case  for  triple-deck  computations.  Sykes 
considered  incompressible  flow  past  this 
shape  as  an  example  of  topographic  effects 
on  the  Earth's  boundary  layer;  nls  results 

are  compared  here  with  results  obtained 

22 

spectrally.  Ragab  and  .'iayfeh  also 
obtained  solutions  for  this  shape  wltn 
moderate  amplitude  (a<2.4),  but  reported 
that  their  finite-difference  program  would 
not  converge  for  a  ■  3 ,  even  with  under- 
relaxation  factors  as  small  as  0.1.  A 

similar  failure  was  reported  for 
22 

.Mapolitano 's  finite-difference  scheme. 

The  spectral  method  described  here  treated 
the  problem  with  no  difficulty,  converging 
in  30  iterations  without  needing  under¬ 
relaxation.  The  spectral  results  are 
shown  in  Figure  5  by  the  solid  curve,  and 
Sykes  results  are  indicated  by  tne  symbols. 
The  grid  parameters  are  N  *  256 , 

AX  •  0.125,  Ny  *  31,  AY  -  0.5  for  the 
spectral  results,  and  Nx  •  256,  AX  *  0.08, 
N y  ■  60,  AY  •  0.25,  for  the  finite-dif¬ 
ference  results  of  Sykes.  The  two  numer¬ 
ical  solutions  compare  fairly  well,  with 
the  spectral  results  displaying  a  smoother 
variation  in  the  reversed  flow  region. 

There  appears  to  be  a  slight  oscillation 
in  the  finite-difference  results,  suggest¬ 
ing  the  need  for  a  finer  grid,  already 
finer  than  that  used  in  the  spectral  com¬ 
putations.  It  is  interesting  that  both 
results  3how  a  sharp  break  in  the  slope 
of  the  snear-stress  curve  Just  downstream 
of  the  leading  edge  of  the  hump.  These 
results  for  surface  stress  in  incompres¬ 
sible  flow  appear  qualitatively  quite 
similar  to  those  for  supersonic  flow  past 
the  ( X )  hump.  The  main  difference  is 
that  the  supersonic  flow  showed  a  stronger 
tendency  to  separate  upstream  of  tne  hump, 
while  the  incompressible  flow  generates  a 
rather  large  separation  bubble  in  the  lee 
of  the  a  »  -3  hump.  These  trends  would 
not  appear  to  contradict  physical  insight 
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regarding  tne  behavior  of  supersonic  and 
incompressible  flows. 

V.  Concluding  Remarks 

This  study  has  demonstrated  that  the 
spectral  method  is  fast  and  effective  for 
solving  triple-deck  problems,  with  the 
ability  to  handle  easily  problems  that 
some  recent  finite-difference  methods  are 
unable  to  solve.  Moreover,  the  accuracy 
of  the  spectral  method  is  compatible  with 
that  of  the  better  finite-difference 
schemes.  In  this  regard,  it  is  often 
stated  tnat  the  spectral  method  is  of 
infinite-order  accuracy,  converging  faster 
than  any  power  of  N .  This  statement  would 
be  true  for  finite  domains  if  the  boundary 
conditions  did  not  Introduce  boundary  dis¬ 
continuities  into  the  function  represented 
by  the  Fourier  series.  Chebyshev  polyno¬ 
mial  representation  is  preferred  for  this 
lb 

reason.  For  infinite  domains,  as  in  the 
present  application,  the  Fourier  Integrals 
themselves  lead  to  algebraic-order  accu¬ 
racy  when  discretized  (2nd  order  when 
using  tne  trapezoidal  rule).  In  addition, 
truncating  the  infinite  range  of  integra¬ 
tion  lead3  to  additional  error  depending 
on  the  asymptotic  behavior  of  the  physical 
and  spectral  functions  Involved.  The  use 
of  too  restricted  a  range  of  integration 
can  lead  to  quite  large  errors.  As  can  be 
seen  from  Eq.  (19),  the  finite-range  (NAX) 
error  in  physical  variables  is  equivalent 
to  grid-size  Aui  error  in  spectral  vari¬ 
ables,  and  vice  versa.  Consequently,  it 
is  necessary  to  use  care  in  choosing  grid 
size  in  both  physical  and  spectral  vari¬ 
ables  to  ensure  an  accurate  solution.  Of 
course,  similar  considerations  hold  for 
finite-difference  solutions. 
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Figure  2.  Surface  Stress  for  Supersonic 
Flow  Past  Contour  H,(X). 
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Figure  3.  Effect  of  Grid  Size  on  Accuracy; 
Supersonic  Flow  Past  Contour  K^fX). 
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Figure  4.  Comparison  of  Supersonic  and 
Incompressible  Flows  Past  Contour  H^(X). 


Figure  5.  Surface  Stress  for  Incompres¬ 
sible  Flow  Past  Contour  H,(x). 


